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[000 1] CHANNEL ESTIMATION METHOD AND SYSTEM USING 

FAST FOURIER TRANSFORMS 

P 

[0002] CROSS REFERENCE TO RELATED APPLICATION(S) 

[0003] This application claims priority from U.S. Patent Application No. 

60/460,852, filed April 4, 2003, which is incorporated by reference as if fully set forth. 

[0004] FIELD OF INVENTION 

[0005] This invention generally relates to channel estimation in wireless 
communications. In particular, the invention relates to low cost channel estimation 
using fast Fourier transform (FFT). 

[0006] BACKGROUND 

[0007] In code division multiple access (CDMA) communication systems, multiple 
communications may be simultaneously sent over a shared frequency spectrum. Each 
communication is distinguished by the code used to transmit the communication. 
[0008] In one type of CDMA communication system, the shared spectrum is time 
divided into frames having a predetermined number of time slots, such as fifteen time 
slots. This system is referred to as a hybrid CDMA/time division multiple access 
(TDMA) communication systems. In another type of CDMA system, uplink and 
downlink communications are restricted to particular time slots. This system is 
referred to as a time division duplex (TJDD) communication system. 
[0009] In a typical TDD/CDMA communication system, communication data is 
sent using communication bursts. Figure 1 illustrates a communication burst 16 
having a midamble 20, a guard period 18 and two data fields 22, 24. The data fields 
22, 24 carry the data of the communication burst 16. The guard period 18 separates 
the communication bursts 16 to allow for a difference in arrival times of bursts 
transmitted from different transmitters. The midamble 20 separates the two data 
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fields 22, 24 and has a known training sequence used to estimate the channel that the 
communication burst 16 experiences. Using the estimated channel response, data from 
the data fields 22, 24 is recovered at a receiver. For the third generation partnership 
(3GPP) wideband CDMA (W-CDMA), based on the burst type, the basic midamble 
codes used to generate midamble bursts have differing lengths. To illustrate, the basic 
midamble code for burst type I has 456 chips while burst type II has 192 chips. 
[0010] FFT is a powerful tool for efficient implementation of a wireless 
communications receiver. One drawback with FFT implementations is that they are 
limited to operating on fields with a predetermined length. 

[0011] It is desirable to have a channel estimator using an FFT engine that has 
the capabilities to handle channel estimation for various types of bursts. 

[0012] SUMMARY OF THE INVENTION 

[0013] In a wireless communication system, channel estimation is performed by 
receiving reference signals having different lengths, processing the reference signals 
using a fast Fourier transform (FFT), and extending the FFT to a desired length L for 
more efficient computation. The FFT is extended to the length L to process a plurality 
of different burst types associated with the reference signals. 

[0014] BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] A more detailed understanding of the invention may be had from the 

following description of a preferred example, given by way of example and to be 

understood in conjunction with the accompanying drawing wherein: 

[0016] Figure 1 is an illustration of a communication burst. 

[0017] Figure 2 is a simplified diagram of a transmitter and a receiver using 

channel estimation. 

[0018] Figure 3 is a system block diagram illustrating the process for 
implementing channel estimation in accordance with a preferred embodiment of the 
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present invention. 

[0019] Figure 4 is a flow chart illustrating the process to calculate F(r) and F(m) 
by extended FFT and divide them element-to-element. 

[0020] Figure 5 is a flow chart illustrating the process to compute F(#*) by the 
extended FFT and then conjugate and scale the result. 

[0021] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 
[0022] Although the preferred embodiments are described in conjunction with a 
preferred TDD/CDMA or TDMA/CDMA communication system, some aspects are also 
applicable to CDMA systems in general, including TD-SCDMA. However, the 
invention in its broad form is envisaged to be applicable to other systems of 
transmission also, without limitation. 

[0023] Hereafter, a wireless transmit/receive unit (WTRU) includes but is not 
limited to a user equipment, mobile station, fixed or mobile subscriber unit, pager, or 
any other type of device capable of operating in a wireless environment. When referred 
to hereafter, a base station includes but is not limited to a base station, Node-B, site 
controller, access point or other interfacing device in a wireless environment. 
[0024] Figure 2 illustrates an embodiment of channel estimation as used in a 
wireless communication system operating in accordance with the present invention. A 
transmitter 30 and a receiver 32 communicate with each other via a wireless radio air 
interface 38. The transmitter 30 may be located at a WTRU or at a base station. The 
receiver 32 may be located at the WTRU and/or the base station. 
[0025] Data symbols to be transmitted to the receiver 32 are processed by a 
modulation and spreading device 34 at the transmitter 30. The modulation and 
spreading device 34 spreads the data with the codes and at a spreading factor(s) 
assigned to the communication(s) carrying the data. The communication(s) are 
radiated by an antenna 36 or antenna array of the transmitter 30 through the wireless 
radio interface 38. 
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[0026] At the receiver 32, the communication(s), possibly along with other 
transmitters' communications, are received at an antenna 40 or antenna array of the 
transmitter 30. The received signal is sampled by a sampling device 42, such as at the 
chip rate or at a multiple of the chip rate, to produce a received vector. The received 
vector is processed by a channel estimation device 46 to estimate the channel impulse 
responses for the received communications. The channel estimation device 46 uses a 
training sequence in the received communication to estimate the channel experienced 
by each communication. A data detection device 44, such as a joint detection device, 
uses the code(s) of the received communication(s) and the estimated impulse 
response(s) to estimate soft symbols of the spread data. 

[0027] The following is a description of an exemplary process for channel 
estimation with a single length FFT. 

[0028] For channel estimation, the received signal r can be expressed as the 
circular convolution of two sequences, the midamble sequence, m and the channel 
impulse response, h by 

[0029] ■ r = m®h Equation (1) 

where ® is defined as the circular convolution operator. 

[0030] In frequency domain, the circular convolution of two signals becomes a 
product of frequency responses of two signals, and the frequency response of the 
resulting output signal becomes 

[003 1] K = K K Equation (2) 

where R is the FFT of time domain signal r , M_ is the FFT of midamble sequence m 
and # is the FFT of channel impulse response h . In short they are expressed by 
R - F(r) , M_ = F (m ) and H_ - F(h) where F( ) is defined as the operator for forward 

FFT and F~* ( ) is defined as the operator for inverse FFT. The "underscore" indicates 
that the signal is a vector. 

[0032] To obtain the estimated channel response, first H_ is calculated by 
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[0033] K = — Equation (3) 

M 

where RjM_ is the element-to-element division of the corresponding two FFT 
sequences. The channel impulse response can be estimated by inverse FFT of II by 
[0034] h = F~ l (H) Equation (4) 

[0035] The whole process can be summarized in one single equation by 
[0036] h = F'\F(r)/F(m)) Equation (5) 

[0037] F(r)/ F(m) in Equation (5) denotes the element-to-element division of FFT 
sequences F(r) and F(m). Note that the forward and inverse FFT are exchangeable in 
the following form as illustrated for vector x: 

[0038] F~ l (x) = j(F(x))* Equation (6) 

where P is the length of FFT. 

[0039] To enable a single-length FFT, the FFT operations used for channel 
estimation as shown in Equation (5) should be extended to a longer desired length. In 
the preferred embodiment, the length of extended FFT should be at least long enough 
to account for all burst types including burst types 1, 2 and 3 so that FFTs in Equation 
(5) will not depend on the burst type, although in other systems and implementations 
the length may vary. To extend a FFT to any longer proper length L, a chirp transform 
algorithm (CTA) is used to compute F(r) and F(m) in Equation (5) by an extended 

FFT. The inverse FFT (IFFT) for F~ l (H) in Equation (4) can be implemented by the 

forward FFT. This can be done by performing FFT on the conjugate signal of # and 

then taking conjugate on the resulting FFT output and scaling it properly as shown in 

Equation (6). CTA is used to compute F(H*) by extended FFT. 

[0040] The procedure for extended FFT and efficient channel estimation is 
described as follows: 

[0041] There are two stages to be performed. The first stage calculates F(r) and 
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F(m) by extended FFT and divides them element-to-element. The second stage 
computes F(H_*) by the extended FFT and then conjugates and scales the result. For 

the following discussions, let P denote the original length of FFT and L be the extended 
length of FFT using tail zero-padding. The original lengths of FFT for example, are 
P=456 and P=192 for burst types 1/3 and 2 respectively. 

[0042] Figure 3 is a block diagram of a system 300 including means for performing 
stages 1 and 2. For stage 1, an element-to-element multipliers 305A, 305B multiply 

the sequences m and r by chirp waveform W"^ for n = 0, 1, 2,...,P-1 where P=456 for 

_ .2tt 

burst types 1/3 or P=192 for burst type 2 and W = e p . In this context, chirp 

waveform is referred to as the waveform generated by W 2 , n = 0, 1, 2,...,P-1. The 
resulting sequences are denoted as m w and r w respectively. A chirp sequence v is 

created such that v = W' { " ™^ for n = 0, 1, 2,. . .,2P-2. Chirp sequence is referred to the 

(n-P+l) 2 

modified sequence based on a chirp waveform generated by W 2 , n = 0, 1, 2,. . .,2P- 
2. Chirp sequence is different from chirp waveform with a shift in index n and with a 
longer waveform length. The chirp transform algorithm refers to the entire process 
that performs the original FFT with extended FFT, in which the chirp waveform and 
chirp sequence are used to transform signals in a proper format suitable for processing. 
[0043] The sequences m w , r w and v are processed by zero padding 310A, 310B, 
3 10C in the tail until the length of the sequences achieves L. Denote the resulting 
sequences as m w z , r wz and v z . L-point FFTs 315,A, 315B, 315C are implemented on 

m w z , r w z and v z each such that F(m w z ), F(r w z ) and F(v z ). Element-to-element 
multipliers 320A, 320B multiply the FFT of m w z and r w z each with FFT of v z such 
that the products are F(m wz ) -F(v z ) and F(r w z )- F(v z ) respectively. L-point inverse 
FFTs (IFFTs) 325A, 325B are implemented on the outputs of multipliers 320A, 320B 
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such that F l (F(m wz )-F(v z )) and F l (F(r wz )>F(v z )) respectively. An element-to- 
element divider 330 divides the outputs of L-point IFFTs 325A, 325B and denotes the 

F~ l (F(r )-F(v ) 

result as H 335 such that H = / y ~ w ' zJ — . Note that only the first P elements 

_ F-\F(m WtZ )F(v z ) 

of sequence H_ are computed and used. 

[0044] Still referring to Figure 3, for stage 2, conjugate device 340 conjugates the 
sequence H that was obtained in the final step of stage 1. Element-to-element 

multiplier 345 multiplies the conjugate sequence H* by chirp waveform W n/ * for n =0, 

1, 2,...,P-1. The result is denoted as H^ w . Zero padding 350 zero pads the conjugate 

sequences H_* w in the tail until the length of the sequence achieves length L. The 

resulting sequence is denoted as h" wz . An L-point FFT 355 is performed on H* wz . 

Element-to-element multiplier 360 multiplies the FFT of H* wz by FFT of zero-padded 

chirp sequence v z such that the product is F(H* W z ) • F(v z ) . An L-point IFFT 365 is 

implemented on the output of multiplier 360 resulting F(H* W z ) • F(v z ) . Element-to- 

* n 2 / 

element multiplier 370 multiplies the sequence F(H_ W Z ) ■ F{v z ) by chirp waveform W h 

for n = 0, 1, 2,. . .,P-1. Note that only the first P elements are calculated and used. The 
output of multiplier 370 is conjugated by conjugate device 375 and the result is scaled 

by scaling device 380 by factor to obtain the estimated channel response. 

[0045] Referring to Figure 4, the procedure for stage 1 by the extended FFT is 
described as follows: 

[0046] In step 405, element-to-element, multiply the sequences m and r by chirp 

waveform W"^ for n = 0, 1, 2,. . . ,P-1 where P=456 for burst types 1/3 or P=192 for burst 
type 2. Denote the resulting sequences as m w and r w respectively. 

[0047] In step 410, create a chirp sequence v such that v = W~ A for n = 0, 1, 
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2,...,2P-2. 

[0048] In step 415, zero pad the sequences m w , r w and v in the tail until the 
length of the sequences achieves L. Denote the resulting sequences as m wz , r wz and 

[0049] In step 420, perform L-point FFT on m wz , r w z and v z each such that 
F(m wz ), F (r wz ) and F(v z ). 

[0050] In step 425, element-to-element multiply the FFT of m w z and r wz each 
with FFT of v z such that the products are F(m wz )- F(v z ) and F(r w z )- F(v z ) 
respectively. 

[0051] In step 430, an L-point inverse FFT is performed on F(m w z )- F(v z ) and 
F(r WtZ )-F(v z ) such that F _1 (F(m wz ) • F(v z )) and F~ l (F(r w z )- F(v z )) respectively. 
[0052] In step 435, element-to-element divide F~ l (F(m wz ) • F(v z )) by 

F" 1 (F(r w z ) • F(v z )) and denote the result as H such that H = . v K ~ w > z) — K -=^- . Note 

F (F(m wz )-F(v z ) 

that only the first P elements of sequence // are computed and used. 

[0053] Note that basic midamble code is fixed in a cell and the chirp sequence is 

constant once it is generated, the computation of F~\F(m wz )-F(v z )) can be pre- 

computed once and for all and stored for each cell. 

[0054] Referring to Figure 5, the procedure for stage 2 is to compute F(H_*) by 
extended FFT and conjugate and scale the result is described as follows: 
[0055] In step 505, conjugate the sequence H_ that was obtained in the final step 
of stage 1 (step 435). 

[0056] In step 510, element-to-element multiply the conjugate sequence H_* by 

chirp waveform W h for n = 0, 1, 2,...,P-1. Denote the result as H_ w . 

[0057] In step 515, zero pad the conjugate sequences l£ w in the tail until the 
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length of the sequence achieves L. Denote the resulting sequence as H_ * w z . 
[0058] In step 520, perform L-point FFT on H w z . 

[0059] In step 525, element-to-element multiply the FFT of H wz by FFT of zero- 
padded chirp sequence v z such that the product is F{i£ wz ) • F(v z ) . 
[0060] In step 530, perform L-point inverse FFT on F(H WZ ) • F(v z ) . 
[0061] In step 535, element-to-element multiply the sequence F~ l (F(H* WZ ) • F(v z )) 

by chirp waveform W^ 1 for n = 0, 1, 2,...,P-1. Note that only the first P elements are 
calculated and used. 

[0062] In step 540, conjugate the result in step 535. 

[0063] In step 545, scale the result in step 540 by factor -i to obtain the estimated 
channel response. 

[0064] The above method of extended FFT using simple zero padding for channel 
estimation is very cost efficient in terms of hardware complexity. High computational 
efficiency is also achievable when the extended FFT length is optimized for specific 
computing algorithms such as a prime factor algorithm (PFA) or Radix-2 algorithm. 
[0065] While this invention has been particularly shown and described with 
reference to preferred embodiments, it will be understood by those skilled in the art 
that various changes in form and details may be made therein without departing from 
the scope of the invention described hereinabove. 



